% solid tide unit test

XSTA(1) = 4075578.385D0;
XSTA(2) =  931852.890D0;
XSTA(3) = 4801570.154D0;
XSUN(1) = 137859926952.015D0;
XSUN(2) = 54228127881.4350D0;
XSUN(3) = 23509422341.6960D0;
XMON(1) = -179996231.920342D0;
XMON(2) = -312468450.131567D0;
XMON(3) = -169288918.592160D0;
YR      = 2009;
MONTH   = 4;
DAY     = 13;
FHR     = 0.00D0;

DXTIDE0(1) = 0.7700420357108125891D-01;
DXTIDE0(2) = 0.6304056321824967613D-01;
DXTIDE0(3) = 0.5516568152597246810D-01;


dxtide=dehanttideinel(XSTA,YR,MONTH,DAY,FHR,XSUN,XMON);
err = dxtide-DXTIDE0

%%%%
XSTA(1) =  1112189.660D0;
XSTA(2) = -4842955.026D0;
XSTA(3) =  3985352.284D0;
XSUN(1) = -54537460436.2357D0;
XSUN(2) =  130244288385.279D0;
XSUN(3) =  56463429031.5996D0;
XMON(1) =  300396716.912D0;
XMON(2) =  243238281.451D0;
XMON(3) =  120548075.939D0;
YR      = 2012;
MONTH   = 7;
DAY     = 13;
FHR     = 0.00D0;

DXTIDE0(1) = -0.2036831479592075833D-01;
DXTIDE0(2) =  0.5658254776225972449D-01;
DXTIDE0(3) = -0.7597679676871742227D-01;

dxtide=dehanttideinel(XSTA,YR,MONTH,DAY,FHR,XSUN,XMON);
err = dxtide-DXTIDE0

XSTA(1) =  1112200.5696D0;
XSTA(2) = -4842957.8511D0;
XSTA(3) =  3985345.9122D0;
XSUN(1) =  100210282451.6279D0;
XSUN(2) =  103055630398.3160D0;
XSUN(3) =  56855096480.4475D0;
XMON(1) =  369817604.4348D0;
XMON(2) =  1897917.5258D0;
XMON(3) =  120804980.8284D0;
YR      = 2015;
MONTH   = 7;
DAY     = 15;
FHR     = 0.00D0 ;

DXTIDE0(1) = .00509570869172363845D0;
DXTIDE0(2) = .0828663025983528700D0;
DXTIDE0(3) = -.0636634925404189617D0;
dxtide=dehanttideinel(XSTA,YR,MONTH,DAY,FHR,XSUN,XMON);
err = dxtide-DXTIDE0

XSTA(1) = 1112152.8166D0; 
XSTA(2) = -4842857.5435D0; 
XSTA(3) = 3985496.1783D0; 
XSUN(1) = 8382471154.1312895D0; 
XSUN(2) = 10512408445.356153D0; 
XSUN(3) = -5360583240.3763866D0; 
XMON(1) = 380934092.93550891D0; 
XMON(2) = 2871428.1904491195D0; 
XMON(3) = 79015680.553570181D0; 
YR      = 2017;
MONTH   = 1;
DAY     = 15;
FHR     = 0.00D0; 

DXTIDE0(1) =  .0050957086917236384D0;
DXTIDE0(2) =  .082866302598352870D0;
DXTIDE0(3) = -.063663492540418962D0 ;

dxtide=dehanttideinel(XSTA,YR,MONTH,DAY,FHR,XSUN,XMON);
err = dxtide-DXTIDE0















